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ABSTRACT 

We extend the confined-density-dependent-mass model to include isospin dependence of the equiv- 
alent quark mass. Within the confined-isospin-density-dependent-mass model, we study the quark 
matter symmetry energy, the stability of strange quark matter, and the properties of quark stars. 
We find that including isospin dependence of the equivalent quark mass can significantly enhance 
the quark matter symmetry energy and stiffen strange quark matter. The recently discovered large 
mass pulsar PSR J1614-2230 with a mass of 1.97 ± O.O4M0 can be well described by a quark star if 
the equivalent quark mass is strongly isospin dependent, indicating that the quark matter symmetry 
energy might be much stronger than the nuclear matter symmetry energy. 
Subject headings: dense matter - equation of state - stars: neutron 



1. INTRODUCTION 

One of fundamental issues in contemporary nuclear 
physics, astrophysics, and cosmology is to investigate 
the properties of strong interaction matter, especially 
its equation of state (EOS), which plays a central role 
in understanding the nuclear structures and reactions, 
many critical issues in astrophysics, and the matter state 
at early universe. Quantum chromodynamics (QCD) is 
believed to be the fundamental theory for the strong 
interaction. Although the perturbative QCD (pQCD) 
has achieved impressive success in describing high en- 
ergy processes, the direct application of QCD to lower 
energy phenomena remains a big challenge in the com- 
munity d ue to the complicated non-per turbative feature 
of QCD ()Fukushima fc Hatsuda I [20lTI ). The ah initio 
lattice QCD (LQCD) numerical Monte Carlo calcula- 
tions provide a solid basis for our knowledge of strong in- 
teraction matter at finite-temperature regime with zero 
baryon density (baryon chemical potential). However, 
the regime of finite baryon density is still inaccessi- 
ble by Monte Carlo becau se of the Fermion sign prob- 
lem ([Barbour et al. 1119860 . 

In terrestrial laboratory, heavy ion collisions (HIC's) 
provide a unique tool to explore the properties of strong 
interaction matter. The experiments of high energy 
HIC's performed (or being performed) in the Relativis- 
tic Heavy Ion Collider (RHIC) at BNL and the Large 
Hadron Collider (LHC) at CERN have revealed many 
interesting features of strong interaction matter at zero 
baryon density and high temperature. Instead of the 
original picture of a hot ideal gas of non-interacting 
deconfincd quarks and gluons at zero baryon density 
and high temperature, the experimental data support 
a new picture that quarks and gluons form a strongly 
interacting system, just like a perfect liquid, in which 
non-p erturbative physics plays an important role ()Tang I 
|2009| ). On the other hand, the properties of strong in- 
teraction matter at higher baryon density regions can 
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be explored by the beam-energy scan program at RHIC 
which aims to give a detailed picture of QCD phase 
structure, especially to locate the so-called QCD criti- 
cal point (jStephanov et al. 1 119981 ). Knowledge of high 
baryon density regions can be further complemented by 
future experiments planned in the Facility for Antipro- 
ton and Ion Research (FAIR) at GSI and the Nuclotron- 
based Ion Cohider Facility (NICA) at JINR. 

In nature, the compact stars provide another way 
to explore the properties of strong interaction matter 
at high baryon density (and low temperature). Neu- 
tron stars (NS's) have been shown to provide natural 
testing grounds of our knowl edge about the EOS of 
neutron-rich nuclear matter (|Lattimer fc Prakash 11200^ 
iSteiner et al. I [20051 ) . In the interior (or core) of NS's, 
there may exist hyperons, meson condensations, and 
even quark matter. Theoretically, NS' s may be con- 
verted to (stran ge) quark s tars ('QS's') (iBonibaci et all 
[200l iStaff elal . 2007; Hcrz og fc Ropke II201H) . which is 
made purely of absolutely stable deconfined m, d, s quark 
matter (with some leptons), i.e., strange quark matter 
(SQM). Although most of observations related to com- 
pact stars can be explained by the conventional NS mod- 
els, the QS hypothesis cannot be conclusively ruled out. 
One important feature of QS's is t hat for a fixed mass, 
QS's have smaller radii than NS's (iKapoor and Shukre I 
1200 If ). It has been argued l[Weber Il2005[ ) that the~mr 
usual small radii exacted from observational data sup- 
port that the compact objects SAX J1808.4C3658, 4U 
1728C34, 4U 1820C30, RX J1856.5C3754 and Her X- 
1 are QS's rather than NS's. The possible existence 
of QS's is one of the most intriguing aspects of mod- 
ern astrophysics and has important implications for as- 
trophysics and the strong interaction physics, especially 
the properties of SQM which ess entially determi ne the 
structure of QS's (llyanenko &: Ku rdgclaidzc 1969: IItoh I 
1970t iBodm cr "1971*: 'Wittcn "1984t iFarhi fc Jafre1[l9M 
Alcock et al. ..1986: Weber 20oC^ 



The EOS of dense quark matter is usually soft due 
to the asymptotic freedom of QCD for quark-quark in- 
teractions at extremely high density. In addition, the 
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EOS of SQM will be further soften due to the addi- 
tion of s quark which contributes a new degree of free- 
dom. Therefore, most of quark matter models pre- 
dict relatively smaller maximum mass of QS's. Re- 
cently, by using the general relativistic Shapiro delay, 
the mass of PSR J1 614-2230 was precisely measured to 
be 1.97 ± 0.04Mo (|Demorest et al.ll20ld) . which gives 
the new record of the maximum mass of pulsars. This 
high mass seems to rule out conventional QS models 
(whose EOS's are soft), although some other models of 
pulsar-like stars wi th quark matter can sti l l describe the 
large mass pulsar d Alford fc Reddv 1 120031: IBaldo et al. I 
200aiRuster fc Rischke „2004HAlford et al. II2005L 12007 
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Bonanno fc Sedra kian 1 120121)" All these models seem 



to indicate that to obtain a large mass (about 2M0) 
pulsar-like star with quark matter, the interaction be- 
tween quarks should be strong, remarkably consistent 
with the finding in high energy HIC's that quarks and 
gluons form a strongly interacting system. 

In QS's, the u-d quark asymmetry (isospin asym- 
metry) could be large, and thus the isovector proper- 
ties of SQM may play an important role. Further- 
more, the quark matter formed in high energy HIC's 
at RHIC/LHC (and future FAIR/NIC A) generally also 
has unequal u and d {u and d) quark numbers, i.e., it 
is isospin asymmetric. In recent years, some interest- 
ing features of QCD phase diagram at finite isospin have 
been revealed b ased on LQCD an d som e phenomeno- 
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the isovector properties of quark matter, which is poorly 
known, especially at finite baryon density. Therefore, it 
is of great interest and critical importance to explore the 
isovector properties of quark matter, which is useful for 
understanding the properties of QS's, the isospin depen- 
dence of hadron-quark phase transition and QCD phase 
diagram, and the isospin effects of partonic dynamics in 
high energy HIC's. 

QS's provide excellent astrophysical laboratory to ex- 
plore the isovector properties of quark matter at high 
baryon density. In the present work, by extending the 
confined-density-dependent-mass (CDDM) model to in- 
clude isospin dependence of the equivalent quark mass, 
we investigate the quark matter symmetry energy and 
the properties of QS's. We find that, although the max- 
imum mass of QS's within the original CDDM model 
is significantly smaller than 2Mq, a strong quark mat- 
ter symmetry energy introduced in the extended CDDM 
model can stiffen SQM and thus the large mass pulsar 
PSR J1614-2230 with a mass of 1.97 ± 0.04Mq can be 
well described by a QS. 

2. THE THEORETICAL FORMULISM 

2.1. The confined isospin and density dependent mass 

model 

According to the Bodrner-\Y itten-Terazawa hypothe- 
sis (iWitten 111983 iWeber II2005I) . SQM might be the true 
ground state of QCD matter (i.e., the strong interaction 



matter) and is absolutely stable. Furthermore, Farhi and 
Jaffe found that SQM is stable near nuc lear saturation 
densi ty for large model parameter space (jFarhi fc Jaffe I 
119841 ). The properties of SQM cannot be calculated 
directly from pQCD or LQCD, because SQM has fi- 
nite baryon density and its energy scale is not very 
high. To understand the properties of SQM, people have 
builded some QCD-inspired effecti ve phenomenological 
models, such as the MIT bag model (jChodos et al. II1974I: 
lAlford et al.l [2005t I Weber I [20051). the Nambu-Jona- 
Lasin i o (NJL) model (IRehberget al.l 119961: IHan et all 



[200lt IRuster fc Rischke I I2004i IMenezes et al. IJ2006D 
the pQCD approach (iFreedman fc Mclerran 1119771 119781: 
iKurkela et al. I I2010D. and the D y son- Schwinger ap- 
proach (iRoberts fc Williams I 119941 : iZong et al. 1 120051: 
IQin et al. 1 120111: ILi et al.l I2011all At extremely high 
baryon densi ty, SQM could be in colo r-flavor-locked 
(CFL) state ([Raiagopal fc Wilczekl [2000t) in which the 
current masses of u, d and s quarks become less impor- 
tant compared with their chemical potentials and the 
quarks have equal fractions with the lepton number den- 
sity being zero according to charge neutrality. 

In quark matter models, one of most important 
things is t o treat quark confi n ement . The MIT 
bag model (iFarhi fc Jaffe I [1981 lAlcock et al.1 [19861) 
and its density dependent versions provide a popu- 
lar way to treat quark confinement. Another pop- 
ular way to treat quark confinement is to vary the 
interact ion part of quark mass, such as the CDDM 
model (iFowler et al I 119811: IChakrabartv et all 1 1989 : 
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Wu et al. ll2008l:lYin fc S u 2008) and th e quasi-particle 
model (ISchertler et al. 1119 97. 1998: Peshier et al. 1120001: 
IHorvath fc LugonesI 120041 : lAlfo^g^eTaL 1 120071) . In the 
present work, we focus on the CDDM model in which the 
quark confinement is modeled by the density dependence 
of the interaction part of quark mass, i.e., the density de- 
pendent equivalent quark mass. 

In the CDDM model, the (equivalent) quark mass in 
quark matter with baryon density tis is usually param- 
eterized as 

D 

ruq = mqo + mi ^ m^o H r, (1) 

ub 

where m^o is the quark current mass and m/ = -^^^ re- 
flects the quark interactions in quark matter which is as- 
sumed to be density dependent, z is the quark mass scal- 
ing parameter, and D is a parameter determined by sta- 
bility arguments of SQM. In the original CDDM model 
used to study two-flavor u-d quark matter (jFowler et al. I 
|1981|) . an inversely linear quark mass scaling, i.e., z = 1 
was assumed and the parameter D was taken to be 
3 times the famous MIT bag constant. The CDDM 
model was later extended to include s quarks to inves- 
tigate the properties of S QM (Chakrabartv et al. "1989t 
lChakrabart7lll99lL [T993. 1996; Benvenuto et al 1995j) . 
Obviously, the CDDM model satisfies two basic features 
of QCD, i.e., the asymptotic freedom and quark confine- 
ment through density dependence of the equivalent quark 



mass, I.e. 
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and lim„ 
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00. For 



two-flavor u-d quark matter, the chiral symmetry is re- 



stored at high density due to linins^aonT-q = if the 
current masses of u and d quarks are neglected. 

The density dependence of the interaction part of the 
quark mass, i.e., mj ~ is phcnomenological in the 
CDDM model, and in principle it should be determined 
by non-perturbative QCD calculations. Instead of the in- 
versely linear density dependence for mj which is based 
on the bag model argument, a quark mass scaling pa- 
rameter of z = 1/3 was derived based on th e in-medium 
chiral condensates and linear confinement (jPeng et al. I 
Il999( ) and has been widely used fo r exploring the prop 
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Zheng et al 1 12004': "P eM et al. 1 12006 1: 1 Wen et al. 
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^^008; Li et al M2011bD . In a recent 
work ()Li et al. II2010I) . Li A. et al. investigated the sta- 
bility of SQM and the properties of the corresponding 
QS's for a wide range of quark mass scalings. Their re- 
sults indicate that the resulting maximum mass always 
lies between 1.5Mq and 1.8Mq for all the scalings cho- 
sen there. This implies that the large mass pulsar PSR 
J1614-2230 with a mass of 1.97 ± 0.04Mq, cannot be a 
QS within the CDDM model. In particular, the maxi- 
mum mass with scaling parameter z = 1/3 is only about 
1.65A/q, significantly less than 1.97± 0.04Mo. 

Physically, the quark-quark effective interaction in 
quark matter should be isospin dependent. Based 
on chir al perturbation theor y, it has been shown re- 
cently (jKaiser fc Weise II2009D that the in- medium den- 
sity dependent chiral condensates are significantly de- 
pendent on the isospin. The isospin dependence of the 
in-medium chiral condensates can also be seen from the 
QCD sum rules (jDrukarev et al. I 120041 : iJeong fc Lee I 
|2012| ) . In addition, the quark-quark interaction in quark 
matter will be screened due to pair creation and in- 
frared divergence a nd the (Debye) sc reening length is also 
isospin dependent (|Dev et al. 1119981 ). These features im- 
ply that the equivalent quark mass in Eq. ([T|) should 
be isospin dependent which is neglected in the CDDM 
model. However, the detailed form of isospin dependence 
of the equivalent quark mass is unknown, and in principle 
it should be determined by non-perturbative QCD calcu- 
lations. In the present work, we extend the CDDM model 
to include the isospin dependence of the quark-quark ef- 
fective interactions by assuming phenomenologically the 
following form for the equivalent quark mass in isospin 
asymmetric quark matter with isospin asymmetry S, i.e.. 



m„ = TO^o + mi + niiso 
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where D/, a, and /? are parameters determining isospin 
dependence of the quark-quark effective interactions in 
quark matter, Tg is the isospin quantum number of 
quarks, and here we set = 1 for g = u (u quarks), 
Tq = —1 for q = d {d quarks), and = for g = s (s 
quarks). The isospin asymmetry is defined as 



(5 = 3 



nd - n.u 



rid 



n,. 
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which equals to ~n:i/nB with the isospin density — 
riu — nd and ub ~ {riu + nd)/'i for two-flavor u-d quark 
matter. The above definition of 5 for qu ark matter has 
been extensively used in the literature (|Di Toro et al. I 



20061 : iPaghara fc Schaffner-Biehch IIMotlDi Toro et al. I 
2OI0I: iShao et al. Il2012[ ). We note that one has (5 = 1 
(— 1) for quark matter converted by pure neutron (pro- 
ton) matter according to the nucleon constituent quark 
structure, consistent with the conventional definition for 
nuclear matter, i.e., ^"^^^ = —nz/riB- In Eq. ([2]), 
the parameters a and /3 should be positive so that we 
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spect the asymptotic freedom and quark confinement. 
In addition, in this work, we set z = 1/3 since an be 
derived based on t he in-medium chir al condensates and 
linear confinement (jPeng et al. Ill999[ ). Obviously, in the 
confined isospin and density dependent mass (CIDDM) 
model, the equivalent quark mass in Eq. ^ satisfies 
the exchange symmetry between u and d quarks which is 
required by isospin symmetry of the strong interaction. 
Therefore, the form of isospin dependence of the equiva- 
lent quark mass in Eq. ([2]) is quite general and respects 
the basic features of QCD. 

2.2. The quark matter symmetry energy 

Si milarly to the cas e of nuclear mat- 
ter (|See. e.g.. Li et all l2008| ). the EOS of quark 
matter consisting of u, d, and s quarks, defined by its 
binding energy per baryon number, can be expanded in 
isospin asymmetry 5 as 

E{nB,6,ns) = EainB^ris) + Esym{nB,ns)d'^ + 0{S'^), 

(4) 

where Eo{nB,ns) = E{nB,5 ~ 0,ns) is the binding en- 
ergy per baryon number in three-flavor u-d-s quark mat- 
ter with equal fraction of u and d quarks; the quark mat- 
ter symmetry energy EsymiriB, ris) is expressed as 



EsymijlB, ris) : 



1 d'^E{nB,S,ns) 



2! 



(5) 
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In Eq. (jU, the absence of odd-order terms in S is due 
to the exchange symmetry between u and d quarks in 
quark matter when one neglects the Coulomb interaction 
among quarks. The higher-order coefficients in S arc usu- 
ally very small. Neglecting the contribution from higher- 
order terms in Eq. (U) leads to the empirical parabolic 
law, i.e., E{nB,5,ns) ~ Eo{nB,ns) + Esym{nB,ris)6'^ for 
the EOS of isospin asymmetric quark matter and the 
quark matter symmetry energy can thus be extracted 
approximately from the following expression 



E. 



-E{nB,6^0,ns)]. 



sym{riB,ns) ~ - [E{nB,S = 3,ns) 



(6) 



In quark matter consisting of u, d, and s quarks, the 
baryon number density is given by ub = {nu + rid + ns)/3 
and the quark number density can be expressed as 



(7) 



where gi = 6 is the degeneracy factor of quarks and vi is 
the Fermi momentum of different quarks (i = it, d, and 
s). Furthermore, the Fermi momenta of u and d quarks 



can be expressed, respectively, as 

Ud = {l + 5/?,)-'y, 



(8) 



where v is the quark Fermi momentum of symmetric u-d 
quark matter at quark number density n = 2nu — 2nd. 
The total energy density of the u-d-s quark matter can 
then be expressed as 
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Using the isospin and density dependent equivalent quark 
masses as in Eq. ([2|) , we can obtain analytically the quark 
matter symmetry energy as 



1 d'^euds/ns 



\5=0 



A + B 
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9m^ 
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(10) 



(11) 



m2 — 3m In 



(12) 



X (L>/7iB"e~'^"«)2, 

and TO = TOiio (or TO.do)+ ^^^^/3 ■ In the present work, we 

assume TOuo = 'Tido = 5.5 MeV and TOso = 80 MeV. In 
the CDDM model, the quark matter symmetry energy is 
reduced to 



EsyminB, ns) 
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It should be noted that the quark matter symmetry en- 
ergy generally depends on the fraction of s quarks in 
quark matter since s quarks contribute to the baryon 
density ub- For two-flavor u-d quark matter, the quark 
matter symmetry energy is reduced to the well-known 

expression, i.e., Esy„,{nB) = ^J+^i ■ 

2.3. Properties of strange quark matter 

For SQM, we assume it is neutrino-free and composed 
of u, d, s quarks and e~ in beta-equilibrium with electric 
charge neutrality. The weak beta-equilibrium condition 
can then be expressed as 



/i„ + fie = = A^s, 



(14) 



where fii {i = u, d, s and e ) is the chemical potential 
of the particles in SQM. Furthermore, the electric charge 



neutrality condition can be written as 
2 1 1 

■^n-u = -Ud + -Us + He 



(15) 



The chemical potential of particles in SQM can be ob- 
tained as 



d£ I — ^ T 5to, driB , I v.. 



driB dui ' \m 
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with 



2.t3 



=V(a;2 + l) + ln (x+x/x^ + l) , (17) 



and e is the total energy density of SQM. One can see 
clearly from Eq. (fT6|) that the chemical potential of 
quarks in SQM has two additional parts compared with 
the case of free Fermi gas, due to the density and isospin 
dependence of the quark mass, respectively. In partic- 
ular, the u quark chemical potential can be expressed 
analytically as 



j—u,d.s 
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TO,,, / \mdj 



[riu +ndY' 
For d and s quarks, we have, respectively. 
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and 



j—u.d,s ■' 
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For electrons, the chemical potential can be expressed as 

(21) 



Me =V 37r2l/e^ + me2. 
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The pressure of SQM can be given by 

j—u^d.s,e 

drrij duB ^ f Vj 
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where —flo is the free-particle contribution and fio can 
be expressed analytically as 



fin 
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487r2 



iy,y^iAi^+Tri^{2vf - 3m^) 



-|-3m^arcsinh 
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Because of the additional parts in the quark chemical 
potentials due to the density and isospin dependence of 
the quark mass, the pressure also has corresponding ad- 
ditional terms. Including such terms is important for 
guaranteeing the thermodynamic self-consistency of the 
model and the Hugcnholtz-Van Hove theorem is then ful- 
filled (Peng et al. 1999). 

3. RESULTS 
3.1. The quark matter symmetry energy 

In Fig. [U we show the baryon number density de- 
pendence of the quark matter symmetry energy in the 
CIDDM model with three typical parameter sets, i.e., DI- 
0, DI-300, and DI-2500. We have considered two typical 
cases, i.e., two- flavor u-d quark matter with ng = and 
u-d-s quark matter with rig = ng. The latter roughly 
corresponds to the situation inside QS's where s quarks 
may have equal fraction as u and d quarks. For the pa- 
rameter set DI-0, wc have Dj = 0, a = 0, /3 = 0, and 
D^/^ = 156 MeV, which corresponds to a typical pa- 
rameter set in the CDDM model. For the parameter set 
DI-300, we have Dj = 300 MeV-fm^", a = 1, /3 = 0.1 
fm^ and D^/^ = 151 McV while Di = 2500 McV-fm^", 
a = 0.8, = 0.1 fm^, and D^^^ = 144 MeV for the 
parameter set DI-2500. 

It is seen from Fig. [T] that the three parameter sets 
DI-0, DI-300, and DI-2500 give very different predictions 
for the density dependent quark matter symmetry en- 
ergy, and thus the three parameter sets allow us to ex- 
plore the quark matter symmetry energy effects. For 
comparison, wc also include in Fig. [1] the nuclear mat- 
ter symmetry energy from the covariant relativ istic mean 
field (RMF) model with interaction NLpS (jLiu et al.l 
I2002D which includes the isovector-scalar S meson field 
and is obtained by fitting the empirical properties of 
asymmetric nuclear matter and describes reasonably well 
the binding energies and charge ra dii of a large num- 
ber of nuclei (jGaitanos et al. II2004D . Although the den- 
sity dependence of the nuclear matter symmetry energy 
is still la rgely uncertainty, especially at supersatur ation 
density ()For a recent review, see, e.g.. Chen Il2012f ). the 
NhpS result nevertheless represents a typical prediction 
of the nuclear matter symmetry energy. One can see 
that the symmetry energy values of two-flavor u-d quark 
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Fig. 1. — (Color online) The quark matter symmetry energy as a 
function of baryon number density in the CIDDM model with three 
parameter sets, i.e., DI-0, DI-300, and DI-2500. The two-flavor u-d 
quark matter with ris = (left window) and the u-d-s quark matter 
with ns = ng (right window) are considered. The nuclear matter 
symmetry energy from the RMF model with interaction NhpS is 
also included for comparison. The symmetry energy values from 
DI-2500 have been divided by a factor of 30 while those of DI-0 
have been multiplied by a factor of 5. 

matter predicted by DI-300 are remarkably in agreement 
with those of nuclear matter with NLpJ while the u-d 
quark matter symmetry energy predicted by DI-2500 are 
about 50 times stronger. On the other hand, the am- 
plitude of the quark matter symmetry energy predicted 
by DI-0 is much smaller than that of the nuclear matter 
symmetry energy. In addition, one can see from Fig. [1] 
that increasing s quark fraction in u-d-s quark matter 
reduces the quark matter symmetry energy as expected 
since s quarks contribute to the baryon density n b while 
the symmetry energy is deflned by per baryon number. 

In the CIDDM model, we can generally increase the 
quark matter symmetry energy by increasing the Di 
value. It should be mentioned that when the value of Di 
parameter is varied, the other three parameters a, /3, and 
D usually need correspondingly readjustment to guaran- 
tee the stability of SQM. As we will see in the following, 
the three parameter sets DI-0, DI-300, and DI-2500 all 
satisfy the stability conditions of SQM, and thus they 
can be used to study the properties of SQM and QS's. 

3.2. The stability of SQM 

Following Farhi and Jaffe (jFarhi fc Jaffelll984[ ). the 
absolute stability of SQM requires that the minimum 
energy per baryon of SQM should be less than the min- 
imum energy per baryon of observed stable nuclei, i.e., 
M(5<^Fe)c2/56 = 930 MeV, and at the same time the 
minimum energy per baryon of the beta-equilibrium two- 
flavor u-d quark matter should be larger than 930 MeV to 
be consistent with the standard nuclear physics. These 
stability conditions usually put very strong constraints 
on the value of the parameters in quark matter models. 

Figure [5] shows the energy per baryon and the corre- 
sponding pressure as a function of the baryon density for 
SQM and two-flavor u-d quark matter in /3-equilibrium 
within the CIDDM model with DI-0, DI-300, and DI- 
2500. One can see that for all the three parameter 
sets DI-0, DI-300, and DI-2500, the minimum energy 
per baryon of the beta-equilibrium two-flavor u-d quark 
matter is larger than 930 MeV while that of SQM is 
less than 930 MeV, satisfying the absolute stable condi- 
tions of SQM. Furthermore, it is seen from Fig. [2] that 



6 




0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 1.5 

He (fm-') 

Fig. 2. — (Color online) Energy per baryon and the corresponding 
pressure as a function of the baryon density for SQM and two-flavor 
u-d quark matter in /3-equilibrium within the CIDDM model with 
DI-0. DI-300, and DI-2500. 
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Fig. 3. — (Color online) Quark fraction as a function of the 
baryon density in SQM within the CIDDM model with DI-0, DI- 
300, and DI-2500 

in all cases, the baryon density at the minimum energy 
per baryon is exactly the zero-pressure density, which is 
consistent with the requirement of thermodynamical self- 
consistency. In particular, we note that the zero-pressure 
density of SQM is 0.24 fm-^, 0.23 frn'^, and 0.21 fm'^ 
for DI-0, DI-300, and DI-2500, respectively, which are not 
so far from the nuclear matter normal density of about 
0.16 fm^^. Moreover, one can see from Fig. [2] that the 
stiffness of SQM increases with the Dj parameter (i.e., 
the quark matter symmetry energy). In addition, we 
have checked the sound speed in the quark matter based 
on the calculated pressure and energy density, and we 
find that the sound speed in all cases is less than the 
speed of light in vacuum, and thus satisfying the casual- 
ity condition. 

In Fig. [21 we show the quark fraction as a function of 
the baryon density in SQM within the CIDDM model 
with DI-0, DI-300, and DI-2500. It is interesting to see 
that the difference among u, d, and s quark fractions 
becomes smaller when the quark matter symmetry en- 
ergy is increased (i.e., from DI-0 to DI-300, and then 
to DI-2500). When the quark matter symmetry energy 
is not so large (i.e., in the cases of DI-0 and DI-300), 
the u, d, and s quark fractions are significantly differ- 
ent, especially at lower baryon densities, which leads 
to a larger isospin asymmetry in SQM. On the other 
hand, a large quark matter symmetry energy (i.e., DI- 
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Fig. 4. — (Color online) Equivalent quark mass as a function of 
the baryon density in SQM within the CIDDM model with DI-0, 
DI-300, and DI-2500. 

2500) significantly reduces the difference of u, d, and 
s quark fractions. In particular, for DI-2500, it is re- 
markable to see that the u, d, and s quark fractions 
become equal and approach a value of about 0.33 for 
Pb ^ 0.4 fm"-^, similar to the results from the picture 
of CFL state. In neutron star matter, the similar sym- 
metry energy effect has been observed, i.e., a strong nu- 
clear matter symmetry energy will give a larger proton 
fraction and thus reduce the difference between neutron 
and pro ton fractions in the beta-eq uilibrium neutron star 
matter (jSee. e.g.. Xu et al. l[2009l ). 

Figure |4] show the equivalent quark mass as a func- 
tion of the baryon density in SQM within the CIDDM 
model with DI-0, DI-300, and DI-2500. It is seen that in 
all cases, the equivalent quark mass increases drastically 
with decreasing baryon density, reflecting the feature of 
quark confinement. Furthermore, interestingly one can 
see a clear isospin splitting of the u and d equivalent 
quark masses in SQM for the parameter sets DI-300 and 
DI-2500; with d quarks having larger equivalent mass 
than u quarks. These features reflect isospin dependence 
of quark-quark effective interactions in isospin asymmet- 
ric quark matter in the present CIDDM model. 

3.3. Quark stars 

Using the EOS's of SQM as shown in Fig. [21 we can 
obtain the mass-radius relation of static QS's by solv- 
ing the Tolman-Oppenheimer-Volkov equation. Shown 
in Fig. [5] is the mass-radius relation for static QS's 
within the CIDDM model with DI-0, DI-300, and DI- 
2500. Indicated by the shaded band in Fig. [5] is the 
latest new holder of the maximum ma ss of pulsars of 
1.97 ± O.O4M0 from PSR J1614-2230 pemorest et all 
l2010f ). For the parameter set DI-2500, we also in- 
clude in Fig. [5] the result for rotating QS's with a spin 
period of 3.15 ms (i.e., the measur ed value for PSR 
J1614-2230), by using th e RNS code (ICook et al \\lM ^ 
iStergioulas fc Friedman I 119951 : iKomatsu et al. I iigs" 
developed by Stergioulas (available as a public dom ain 
program at http : / / www . gravity, phys . uwm . edu/ rns / j ) . 
The radius value of the rotating QS's is taken at the 
equator. 

From Fig. [5l one can see for the parameter set DI-0 
which corresponds to the case of the CDDM model and 
has much small quark matter symmetry, the maximum 
mass of the static QS's is 1.65M0 and the corresponding 
radius is 9.60 km. Therefore, the maximum mass with 
DI-0 is significantly smaller than the observed value of 
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Fig. 5. — (Color online) Mass-radius relation for static quark stars 
within the CIDDM model with Dl-0, Dl-300, and Dl-2500. The 
result for rotating quark starts with a spin period of 3.15 ms is also 
shown for the ease of Dl-2500 with the radius at the equator. The 
shaded band represents the latest new holder of the maximum mass 
of pul sars of 1.97±O.O4M0 from PSR J1614-2230 lDemorest eFall 
120101) . 

1.97± 0.04Mq for PSR J1614-2230. For the parameter 
set DI-300 whose prediction of quark matter symmetry 
energy has almost same ampUtude as the nuclear matter 
symmetry energy, the maximum mass of the static QS's 
is increased to 1.78Mo with the corresponding radius 
of 10.40 km, indicating that increasing the Dj parame- 
ter (and thus the quark matter symmetry energy) in the 
CIDDM model can significantly enhance the maximum 
mass of the static QS's. The Dl-2500 parameter set with 
Dj = 2500 MeV-fm^" gives rise to a maximum mass of 

I. 93Af0 for static QS's with the corresponding radius of 

II. 12 km, which is consistent with the observed mass 
of 1.97 ± 0.04Mq for PSR J1614-2230 within the error 
bar. Furthermore, considering the rotation of QS's with 
a spin period of 3.15 ms as the measured value for PSR 
J1614-2230, we obtain the maximum mass of 1.97Mq for 
rotating QS's with the corresponding radius of 11.33 km 
for the parameter set Dl-2500, which is nicely in agree- 
ment with the observed maximum mass of 1.97±0.04Mq. 
Our results imply that in the CIDDM model, a larger 
Di parameter (and thus a strong isospin dependence of 
the equivalent quark mass) is necessary to describe the 
observed large mass of 1.97 ± O.O4M0 for a QS, which 
means the amplitude of the quark matter symmetry en- 
ergy should be much larger than that of the nuclear mat- 
ter symmetry energy. These features indicate that PSR 
J1614-2230 could be a QS in the CIDDM model, and if 
PSR J1614-2230 were indeed a QS, it can put important 
constraint on the isovector properties of quark matter, 
especially the quark matter symmetry energy. 

For an EOS of SQM, it is physically interesting to 
determine the maximum mass of QS's at the maxi- 
mum rotation frequency constrained by the mass shed- 
ding and the secular instability with respect to axisym- 
mctric perturbations, which essentially corresponds to 
the maximum mass of rotating QS's that the EOS can 
support. The maximum angular frequency flmax of 
a rotating QS can be obtained from its static mass 
Mq °* and ra dius by using t he empirical formula 

proposed by (|Gourgoulhon et aL~lll999|l , i.e., i^max = 



TABLE 1 

The maximum mass, the corresponding radius and central 
baryon number density of the static quark stars, the 
maximum rotational frequency fmax for maximum-mass 
static quark stars as well as the corresponding 

gravitational mass and equatorial radius at fmax- within 

THE CIDDM MODEL WITH Dl-0, Dl-300, AND Dl-2500. 





DI-O 


Dl-300 


Di-2500 


M/MQ[static) 


1.65 


1.78 


1.93 


R{km) (static) 


9.60 


10.40 


11.12 


Central density(fm~^) 


1.31 


1.11 


1.06 


fmax (Hz) 


1680 


1547 


1458 


M/Mq (at fmax) 


1.78 


2.12 


2.43 


R{km) (equator at fmax) 


9.93 


11.6 


14.2 



773O(M|)*°VM0)i/2(i?|Jg/lOkm)-3/2 rad-s"!. 

In Table [1] we list the maximum rotational frequency 
fmax for maximum-mass static QS's as well as the cor- 
responding gravitational mass and equatorial radius at 
fmax, within the CIDDM model with DI-O, DI-300, and 
DI-2500. For completeness, we also include in Table [T] 
the results for the maximum mass, the corresponding 
radius and central baryon number density of the static 
QS's. From Table [TJ one can see the maximum rota- 
tional frequency fmax decreases with Dj while the cor- 
responding mass and equatorial radius increase with Dj . 
In particular, for the parameter set DI-2500, we obtain 
fmax = 1458 Hz, and the corresponding mass is 2.43M0 
with radius of 14.2 km, which essentially corresponds to 
the maximum mass configuration of rotating QS's that 
the parameter set DI-2500 can support. 

4. CONCLUSION 

We extend the confincd-density-dependent-mass model 
in which the quark confinement is modeled by the 
density-dependent quark masses to include isospin de- 
pendence of the quark mass. Within the confined- 
isospin-density-dependent-mass model, we study the 
quark matter symmetry energy, the stability of strange 
quark matter, and the properties of quark stars. We 
find that including isospin dependence of the quark mass 
can significantly change the quark matter symmetry en- 
ergy, and the recently discovered new holder of pulsar 
maximum mass, i.e., the millisecond PSR J1614-2230 of 
1.97 ± O.O4M0, can be well described by a quark star 
with a strong isospin dependence of the quark mass, cor- 
responding to a strong quark matter symmetry energy. 
Our results indicate that if PSR J1614-2230 were indeed 
a quark star, it can put important constraint on the 
isovector properties of quark matter, especially the quark 
matter symmetry energy. In particular, we find that PSR 
J1614-2230 can be well described by a quark star within 
the confined-isospin-density-dependent-mass model with 
the parameter set DI-2500, indicating that the quark 
matter symmetry energy might be much stronger than 
the nuclear matter symmetry energy. 

In the present work, we have mainly focused on 
the quark matter symmetry energy and the proper- 
ties of quark stars within the confined-isospin-density- 
dcpendcnt-mass model. In future, it will be interesting 
to see how the present results change if other quark mat- 
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ter models are used, and how the isovector properties 
of quark matter, especiaUy the quark matter symmetry 
energy, will affect other issues such as the quark-hadron 
phase transition at finite isospin density, the partonic 
dynamics in high energy HIC's induced by neutron-rich 
nuclei, and so on. These works are in progress. 
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